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Abstract 

We study the asymptotic diffusion processes with (generally nonlocal) open boundaries in one 
dimension which are exactly solvable by means of the recently developed recursion formula. We 
investigate the stationary states, which cannot be determined in an elementary way. We give the 
equation which includes an auxiliary parameter and determines possible boundary conditions for the 
model to be solved exactly. With the help of that equation, we analyze the density current and 
concentration. We classify the phases according to them. 

1 Introduction 

The study of the systems with reaction and diffusion in one dimension has been an attracting problem, 
or example, models including both reaction and diffusion have been considered in view of chemical 
processes like A + B — > C + D, where A,B,C and D represent some sorts of molecules. Models with only 
pure diffusion are not trivial, either. They have been studied in relation with interface growth, traffic 
flow and so on|l[. Some experiments corresponding with the recent theoretical developments have been 
proceeded 0]. Recently, a problem concerning protein synthesis has been considered in this scheme Q. 
To investigate reaction-diffusion models, a variety of methods have been used. Once taking a continuum 
limit, traditional technique of the field theory such as renormalization group method, mean-field or 
variational approximation are available ||, [|. 

In this paper, we consider the models on a one-dimensional discretized lattice. Some of the models 
on a lattice are shown to be integrable by means of the Jordan- Wigner transformation, Bethe ansatz(8| 
and so on. 

The definitions and notations of reaction-diffusion models considered in this paper are as follows. We 
consider a one-dimensional lattice, where lattice sites are numbered 1,2,- • -,L{= the number of sites), and 
several sorts of particles living on the lattice. The total number of particles is not larger than that of 
lattice sites, i.e. empty sites represented by <fi are allowed. The stochastic process obeys the following 
continuous flow chart. 

1. At an instant t, a pair of neighboring sites, say the i-th and (i + l)-th sites, is chosen randomly. 
We assume that the sorts of particles on the i-th and (i + l)-th sites are A and B, respectively. 

2. The particles on i and i + 1 (A and B) react into the same or other sorts of particles C and D 
respectively in an infinitesimal time interval dt with a given probability w l c _ A D _ B (C, D)dt. This 
process is represented by A + B ^ C + D and we must tell the difference of the order of particles 
(or sites), i.e. A + B is different from B + A in general. 

3. The time should be set t — > t + dt. 

4. A process is repeated until t amounts to the final time tj defined properly. We often put 
tf = +oo. 
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In this paper, we concentrate in models including only particles of one species, A, and the parameters 
fj,,i/,a and (3 in w l a p(p, v) being Z2 variables. Because of the preservation of the probability, we have the 
relation 

^0,00^)= w l r , s {v-r,v - s), (1) 

r,s 

hereafter means that n\ = ni = ■ ■ ■ = in the summation is excluded. Let us introduce the 

m,ri2,— 

probability distribution Pi (n , • • • , ti ; t) meaning the probability to find the i-th site to be occupied by 
A (empty) if Ti = 1 (0) at the time t. The time evolution equation of Pi is a linear one denoted as 

d 

Q- t P L (Tu---,T L ;t) = (HPi)( Tl ,---,ri;t), (2) 

where H is a linear operator called Hamiltonian on the analogy of the Schrodinger equation. Needless 
to say, H can be written in terms of the local interactions characterized by w's as we will see later. We 
consider the ground state, a solution of HPi = 0, which the system will reach after a long time. 

Especially, we deal with problems of the stationary state of an open chain with pure asymmetric ex- 
clusion processes (ASEP). Although the setting is quite simple, to obtain the exact form of the stationary 
state is not easy. With some classes of models, useful recursion relations to determine physical quantities 
were indicated by Derrida, Domany and Mukamel (DDM)[||]. They solved the recursion relation for a 
special case and predicted the phase diagram by means of the mean-field approximation. The solutions 
with general local boundary condition were obtained by Schiitz and Domany (SD) without using the 
mean-field approximation [{Tcfl . 

In this paper, we will consider general boundary conditions including nonlocal boundary interactions 
so that a recursion relation similar to DDM's might be available, i.e. integrable. We calculate some 
physical quantities, the partition function, density current and concentration (one-point function), exactly 
and classify the phases of the system. The plan of this paper is as follows. In the following §Q, we have 
a brief sketch of the previous works by DDM and SD. An equation for the possible integrable boundary 
conditions is given in §[|. We calculate the partition function, density current and concentration in §[| 
and §||. An unknown example is analyzed in §|| Section is left for the summary and conclusion. 



2 Asymmetric Exclusion Process 

In this section, we consider pure-diffusion systems with boundaries, i.e. open chains. Although it seems 
to be very simple, even determination of the ground state is not easy. An epoch-making work in this 
field was done by DDM[||]. They made use of a recursion relation among the ground state with different 
lattice length, and obtained the solution for a special case (a = (3 = 1 in the later notation). The phase 
diagram classified by the strengths of the boundary injection and subtraction rates were conjectured by 
means of the mean-field approximation. SD proved DDM's conjecture to be true by obtaining the exact 
form of the solutions for the general boundary conditions. Another remarkable succeed is development 
of the so-called matrix product method introduced in ref . [JO]] . This method enables correlation functions 



to be calculated algebraically and exactly. Several extension can be found in refs.[12, [b|, 14, 15, 16], for 
example. Recently, this method has been applied for the models including coagulation and decoagulation 
as well as pure diffusion! 17]. 



In this section, we concentrate in the simplest ASEP model on a lattice with length L |9[ defined as 
follows. 
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Figure 1: Pure diffusion model with open boundaries. 

1. Left edge {% = 1) 

At the left edge, the injection and extraction of a particle are possible 

<E> — ► A with a rate w\(l) = a 

or 

A — > <E> with u>i(0) = 7. 

2. Bulk property (« = 1, • • • , L) 

Only diffusion is allowed in the bulk like 

A + ®^® + A with lOn(Ol) = p 

or 

® + A^A + $ with w n (10) = q. 

3. Right edge (i = L) 

At the right edge, the injection and extraction of a particle are possible as well as at the left edge. 

^ A with w[(l) = 5 

or 

A^§ with ^i(O) = /3. 
We assume the conservation of the probability, for example 

$ — > $ with I - a 

at the left edge. The Hamiltonian of this model can be shown to be equivalent to that of XXZ spin 
chain with integrable boundary terms Jl^, [l^]. In this section, we consider the case with q = 7 = 5 = 
for simplicity. Considering the rescaling of the time, we can put p = 1 without any loss of the generality. 
We are interested in the stationary solution of the equation 

= J2 W^K^r'.T!) 

o-i=0,l 
L-l 

+ Y1 K^T\ + + \ P L{Tl, TJ-i, ffi, 0-i+i, Ti+2, • • • , T L ) 

1 = 1 <Tj,<Tj+l=0,l 

+ J2 ih')^ L P L {r u ---,T L - U a L ), (3) 

a L =0,l 
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where 



hi 



h' 



a —7 \ _ I ot 
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5 -p\ = ( Q -p 

-s a ) \ p 



(4) 



(5) 



(6) 



The readers should remember that we do not have to restrict parameters less than one once we write 
the master equation in the above form by rescaling the time. We can consider that they are arbitrary 
positive numbers. 

In DDM's and the following work, it has been found that that the solution of eq.(^) with general 
p,g,a,/?,7 and 5 satisfies surprising recursion relations for L as follows. 



(Mti^ "!^, ~',tl) = -x Ti P L -i{t 2 , ■ ■ -,tl), 

cri=0,l 
L-l 

8=1 (7i,(7j + l=0,l 

= X T .PL-l(n, ■ ■ ■ ,Tj_i,Tj + i, ■■■ ,T L ) - X n+1 PL-l(Tl, ■ ■ ■ ,Tj,Tj +2 , 

12 ( h ')TL P L(n, ■■■ , tl-i, ■-,<7l) = Xt l Pl-i{t~i, • • • , 7L-l), 

<T L =0,1 



(7) 



(8) 
(9) 



where xq and x\ are some numbers, which can be set as — xo = x\ = 1. It should be remarked that the 
normalization property of the probability, -Pl(°i> " * * j °x) = 1 is no more correct. We must renormalize 
the distribution function after obtaining its form||, |irj]. We consider the normalization, i.e. partition 
function, in §Q in detail. 

With the help of this recursion, we can calculate the correlation functions, for example the concen- 
tration, with arbitrary L in principle. Contrary to the calculation of the concentration, those of general 
n-point function are difficult, in fact. Some technology to calculate the general correlation functions like 
the matrix product method\ll] have been developed. 



3 General Boundary Condition 

As we have seen in the previous section, the recursion relation (0)-(P) will play a crucial role to solve the 
model exactly. We consider this recursion relation in detail. In general, we consider a system denoted 
symbolically by 

n 

, - 

Q — o — o — o — ••• — o — o, n = 0, 1,2,.... 
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Here we introduce an s-valued variable X(= Pi, ... , P s ) living on O a t t ne left edge and a usual 2-valued 
variable a(= 0, 1) on o. Let us denote X £ V s and a G V 2 from now. The connection — o — o — represents 
a pure asymmetric hopping, which we put p = 1 and g = in this section. At the right edge (—0), a 
particle is extracted with (3 = 1 (,6 = 0). 

Denoting the distribution function for the ground state as f n (X, uq, ■ ■ ■ , cr n ), the master equation 
takes the following form 

Yev 3 po 

n—l 

+ E h afaZ\fn(X,a , . . . , <?i-\, pi, p i+X , (T i+2 , . . . ,(T n ) + (h!) P a n J n (X,(TQ, . . .,a n -l,p n ) 
i=0 

(10) 

with the boundary matrix H X<J , which will be determined later. Other matrices h and h are defined in 
eqs.(||) and (||) with p = l,q = 0,(3 = 1 and 6 = 0. 

For example, let us consider the simplest ASEP model to make the notation clear. In this case, we 
can set X £ V2 (i.e. Pi = 1, P2 = 0) and 

Ki = Ki + 

It is obvious that we can reformulate the simplest ASEP model with X G V 2 i. Anyway, that model 
can be considered as the simplest example in the following formalism. We want to answer the following 
question. 

Problem: What is the condition for the recursion relation 

f n (X,po,pt, . . . ,p n -i, 1) = f n -l(X,p ,pi, . . . ,Pn-l) for n = 1,2,... ? (11) 

This relation comes from the recursion relation @. We look for the answer to the above question. To 
begin with, we write down the master equation ( JlOj ) explicitly. For n = 0, we can easily obtain 

]T H Y x {f (Y,p) + f Q (X,l) = 0, (12) 
ygv s ,p=o,i 

£ H x p f (Y,p)-f (X,l) = 0, (13) 
Yev 3 ,p=o,i 

which corresponds to the initial condition for the distribution function in the following consideration. 
Next we consider the equation for n = 1; 

£ H^;if(Y,a 1 ,p 2 ) + h^f(X,a u a 2 ) + (h')%f(X, Pl ,a 2 ) = 0. 

Y,a!_ 

The independent elements of the above equation are the following three, 

Y,H x Ui(Y,o-,0) + f (X,0) = 0, (14) 

Y,a 

J2H%fi(Y,a,0)-f (X,0) = 0, (15) 

Y,a 

h(X,l,0) = f (X,l) + f (X,0), (16) 
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where we made use of the relation (|ll|) and eq.(|l6[) takes the same form as that for the simplest ASEP 
model. Similarly, for general n(> 0), we obtain the equation 

^^/ n (y,a,0,0,...,0) + / n _ 1 (X,0,...,0) = 0, (17) 

Y,a 

Y / H^ f n (Y,a,0,0,...,0)-f n - 1 (X,0,...,0) = (18) 



and 



fn\X, (To, fJj_i, 1, 0, CJi_|_2, CT n ) 

= fn-l\X, (Tq, ... , 1, C"j +2 , . . . , a n ) + f n -i(X, (To, ... , <Tj_i, 0, (Tj + 2, . . . , CT„ 



(19) 

for % = 0, . . . , n — 1. Upon setting new variables 

^(X) = / n (X,0,0,...,0), y n (X) = / n (X,l,0,...,0), n = 0,1,2,..., (20) 
we obtain the following equation 

J2H] c t ixn(Y) + J2 H x 1 iyn(Y) + x n - l (X) = 0, n = l,2,..., (21) 

Y Y 

Y,Hl° Q x n {Y)+Y,Hlly n {Y)-x n ^{X) = 0, n = l,2,..., (22) 
y y 

J2HExo(Y)+Y, H x\yo(Y) + yo(X) = (23) 



y y 



^^o x (y) + 5:^o(n-yoW = (24) 



y n (X) = x n _ip0+j/ n _ip0, n = l,2,..., (25) 



which is the condition for the recursion relation (|TT|) . 

We consider the simplest ASEP model again as an example to make the notations clear. Now we put 
X G V\ instead of V2 and can set H x p a = (hi)?, where (/ii)o = «, (hi)® = —a (hi)\ = 7 and (/ii)o = —7 
at the left edge. We can ignore the argument X in x n and y n in this notation. The equation (^l|)-(p4|) 
reads simply 



a 


7)1 


•En 


—a 






a 


7)' 


f ^0 


—a 




[ yo 



n=l,2. 



Therefore, we obtain 



■Era— 1 
%n— 1 

-yo / 



1 7 

:co = 7 + l> ,2/0 = «, x n = + -y n+ i, 



which is compatible with the relation from (p5|), 

yn — -Era— 1 + 2/ra— 1, 
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therefore, the wave functions in the simplest ASEP model satisfy the relation (|ll]) as already known. 

Next let us eliminate the distribution functions in eqs.(21)-(p5|) and see the condition including only 
the parameters in the boundary matrix H. For this sake, we introduce an auxiliary variable t, functions 



vectors 



v n 



<t>(Pf,t) 



j=0 

oo 

3=0 



*(x n (P 1 ),y n (P 1 ),.. 
*(0(P i; f),V(Pi;t), 



.,0(P s ;f),V(P s ;i)), 



and (2s) x (2s) block diagonalized matrices 

Ei = 



ei 



^0 

/ e 2 



e-2 





ei y 

0^ 



ei 



-1 
1 







Pi 








We want to solve the recursion relation (|2 



d2 j 

|25|). Equation (p5|) reads 



and eqs. (|21| 



mean 



^(P;t) + (l-t)^(P;t) = y (P) 

= -Eiv n , 
Hvq = -E 2 v , 
H£(t) = -tExi(t) - E 2 v . 



(26) 

(27) 
(28) 
(29) 



The initial value yo(Pi) can be determined from ( |28| ) by the determinant of the (2s — 1) x (2s — 1) comatrix 
of H + E 2 like 





H Pl ° 

rrPlO 
H Ptl 


+ i ■ 


u PxO 

ttP,0 

• M P 1 1 


rr-Pi+lO 

H P0 

ttPI+iO 

H P 1 1 


■ H Ps0 

n P!0 

TtPsO 

■ H P 1 1 


ttPsI 
tl P 1 
ttPsI 

H P 1 1 


yo(Pi) = 


ttPiO 
H P,0 

ttPiO 

n P t +iO 


ttPiI 
H P l 

ttPiI 

U P 1+ iO 


ttP 1 
• H PiO 

TTPtO 

' n Pi+iO 


TjPi + lO 
TJ-Pi+lO 

n Pi+lO ' 


TtPsO 

• H P l 

TjPsO 

' U P 1+ iO 


TjPsl 

H P>0 
ttPsI 
n Pi+lO 




rrPlO 
U P S 
ttPiO 
H Psl 


H Pl1 

ttPiI 
H Psl 


n PsO 


H PsO 
H Psl 


■ H p °° 

U PsO 
TjPsO 
• H Psl 


H P pfo - 1 

h p : i + 1 
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Therefore, the condition for the recursion relation is the same as the constraint for the boundary matrix 
that there must be a solution for the simultaneous equations (|2^)-(|29|). Let us introduce a (2s + 1) x 
(2s + 1) matrix Mij(t) for i and j {i, j = 1,2,... , s), 



Mij{t) 



where (2s — 1) x (2s) matrix [H 
(H + tEi) like 

(H + tEx)' ~- 





Vo(Pl) 


(H + tEtY 


-Vo(Pl) 




yo(Ps) 




Vo(Pi) 


\-te Pj0 + (1 - t)e Pjl ) 


yo(Pj) , 



\ 



(30) 



/ 



tE{)' is obtained by eliminating the (P s l)-th column vector from 



( H 
H 



PiO 



1 H P 1 



PiO 

Pll T L -"Pll 



H 
H 



PsO 

PiO 

PsO 

Pll 



H 



Pal 

Pll 



PiO 

PsO 



H 



Pii 

PsO 



H 



PsO 
PsO 



t H 



Psl 
PsO 



and ep.j is the unit vector in the (Pjj)-direction, for example, 



e Pl0 =* (1,0,..., 0), e Pl i =* (0,1,..., 0). 

Needless to say, there exist the other equivalent forms of Mj,-, which are obtained by eliminating (Pjl)-th 
(1 < i < s) column vector from {H + tE\). With the above matrices Mjj(t), the constraint can be written 
in the determinant forms as 

rriijit) = det Mij(t) = (31) 
for all % and j = 1, 2, . . . , s) and arbitrary t. The above equations are the main results in this section. 



4 Calculation of Partition Function 

To calculate physical quantities such as the density currents and concentration, we have to obtain the 
form of the partition function Z n : 

z n= X! fn(X,a Q ,ai, . . . ,a n ). (32) 

X = Pl,...,P s {CT i= 0,l} 

Several quantities can be extracted from this function. We calculate the partition function Z n with 
general n exactly by means of a similar technique to that used in DDM's and SD's papers^, |lQ|[ . First, 
let us assume that we have already solved the master equation for some small n's. It is nothing but a 
diagonalization problem for small dimensional matrices. Moreover, let us accept the following assumption. 

Assumption: With proper constants Sj and ai (i = 1,2, ... ,1), the following recursion 
relation is satisfied. 

l 

Y, f n (X,0,...,0)- ]T / n+a (X,0,...,0,0)=5>i<i?- (33) 

X=Pl,...,P s X=Pl,...,P s 8=1 
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Needless to say, the above relation should be checked in the following context. Although we have not 
succeeded in a general proof, we can show the above relation for specific models inductively. In the 
general treatment proceeded in this section, we provide that this relation is true. We prepare constants 
u and v defined by 



u = £ / (X,0)-£^, 
x=p 1 ,...,p s i a% 

X=Pi,...,P„ 

Let us introduce the following value Y n (k), a similar quantity to HLP's in §^, defined by 

Y n{k) = H (1 - Tn)(l - r n _i) • • • (1 - T fc )/„(X,r ,ri, . . . ,r n ), 

x&Vs { n =o,i} 

n>0,0<k<n, 

Y n (n + l)= Y fn{X,TQ,Ti,...,T n ), n > 0, 

XeVs {n=o,i} 

K_i(0) = u, 
and auxiliary functions L p (X) as 

oo 

L p (X) = Y^ n+1 Y n (n-p), p = -1,0,1,... 

n=p 

with an extra parameter A. By definition, we find 

oo oo 

L-l(A) = ^ n+1 Yn(n + l) = u + Y^Zn. 

n=— 1 n=0 

It is not difficult to prove the following relations 

(1-A)L_!(A)-L (A) = u + vX, (34) 
L P+1 (A) - L P (X) + AL P _!(A) = A ?+1 (Vi(0) " Y p (0)) 

= — 5Z s i a f~\ f° r P > 0. (35) 



Now, we solve the recursion relation (35) under the initial condition (|34l) and the condition that the first 
term in L p (X) is proportional A p+1 . 

For this sake, introducing a generating function w(x, A) with another extra variable x as 



w[x, 



p=0 

the recursion relation ( |35|) becomes 



1 - x + Xx 2 )w{x, A) = L_i + x(L - L_i) + ^ x p+1 A p ^ s;a: 

p=l V i 



p-1 
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giving the solution 



n+2 n+2 n+1 n+1 

u l ~ U 2 t U l ~^2 it r n 

H (+o — L-i 



U)\ — U>2 



UJl — U>2 



uJ 



n+1 



+ 



UJ. 



n+1 



+ 



OiA) 



n+1 



with functions 



(UJl — U)2)(UJ1 — OjA) (UJ2 — U)l){uJ2 — CLi\) (diX — UJl)((liX — U)2) 



1- V1-4A 1 + VI - 4A 

0J\ = , U>2 = • 



Demanding the general condition L p (X) ~ A p+1 , we obtain a simple equation 



-uj 2 L-i - (L - L_i) - ^2 



1 



ll>2 — atX 



and the initial condition ([3). Finally, we see 



i 



+ vX- ^SiX- 



1 



2 / , , \ 2 /, , \ 3 



ti 1 

A 



-A 



S,;a ? ; — A 



A 



3 / 



A 



Expanding L_i(A) as for A with helpful relations 

A/ ~~ n!(n + m)! 

n=0 v ' 



A n , m = l,2, 



gives the expressions 



v n + 1 



2 2n+ 1. 
(2n + 2)! 



Z (o) + Z (i) 



- 

(n + 1)! (n + 2)!' 



i=l 



with a function R n {a) defined by 



Rn{0) = E( m + 2 ) 



m=l 



(2n - m + 1)! 
(n + 2)! (n - m)!* 



(36) 



(37) 



(38) 



The above is the exact form of the partition functions. 

Let us investigate the asymptotic behavior of R a {x) defined in ( |38| ) with large n for later convenience. 
We define functions 

(2m — m + 1)! m 

r n ,m(fl) = (m + 2) 7 — — -to , m=l,2, ...,n. 

(n + 2)! (n — mj! 
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• Region a < 2 

It is not difficult to see that the asymptotic behavior of the function r nm {a) for large n is 

r n ,i{a) = ^aZ^\ r n>m (a) ~ m (^) r n ,i(a). (39) 

Therefore, the asymptotic form of the function R n {a) for large n can be written as 

Rn(a) ~c(o) x 4°), (40) 
where the function c(a) is independent of n 

• Region a > 2 

For large n, because R n {o) is dominated by maxi<k<n ?*n,ft(a), we obtain 

fln(a) ~ -^(a)", (41) 

a — 2 i i a— 2 

fa-2\-*=i / a-2\ 1+ — fAn . 
r(a) = a- (— j (l + — j . (42) 

We will make use of the above expression in the next section. 

5 Density Current and Concentration 

In this section, we calculate some physical quantities and draw some phase diagrams, which just means 
the diagrams of the difference of the asymptotic forms of the density current and one-point function. 
First, we consider the density current 

j _ j- m fn(x, op, • • • , at-i, l? Q; &k+2, • • • , gvQ 

~~ n-^oo 2-^2-^2.^ <7 ' 

X j^k {a} 

which is actually independent of k{= 0, 1, • • • , n) and corresponds to 

j= Hm %±. (43) 

From the asymptotic form of the function R n (a), ( p0[ ) and (f4T[), we can calculate the current as follows. 
It must be remarked that the current does not depend on but only on the eigenvalues a^. 

• Case 1: aj > 2 for some i's. 
Let us assume 

o-l > 0,2 > ■ ■ ■ > o-k > 2 > at+i > ■ ■ ■ > ai. 
From the asymptotic form of R n (a), it is not so difficult to show 

\ai-l) \ ai-lj 

We must remark J < 1/4 in this case. 
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L 

o 



0.0 

i/L 



1 



Figure 2: Concentration ((rj)) in the low density phase, a and f3 correspond to those in the simplest 
ASEP model with p = 1,7 = 5 = 0. This notation is also used in Figs. ||] and |j. 



(n) 0.5 







a + [1 > 1 - 




a + /3 < 1 ■ 



i/L 



Figure 3: Concentration in the high density phase. 



• Case 2: a, < 2 for all £(= 1, 2, . . . , I). 

According to eq. ([!(]), the partition function Z n is proportional to Zn^ in this case. Therefore, we 



can show 



J = lim 



Z 



(o) 

n-l 



,(0) 



1 

4' 



Hence, the density current in this region is not so sensitive to the parameters as that in the previous 
region. We call this region maximal current phase^, [To| . 

We can say nothing more from only the density current. However, the above former region is classified 
further into several phases. Let us consider the cases with two eigenvalues, i.e. / = 2. In the cases with 
two eigenvalues, the concentration, i.e. one-point function, 



(a, 



clef 



1 



n—mn 



Zn 



^ n _ m ,l/n(A,cT ,tTl,...,tT n ) 



X=P 1 ,...,P S {cr i= 0,l} 



z 



1 " m 



(44) 



separates the low and high density phases^], pX| ]. The shapes of the concentrations in the low and high 
density and maximal current phases are shown in Figs.|2],|| and |], respectively. As shown in a general 
scheme done in []l~0| , without the mean-field approximation, systems with just two eigenvalues can 
undergo no other phase transition other than those in Figsj2||4|. The problems on systems with more than 
three eigenvalues should be open. 
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Figure 4: Concentration in the maximal current phase. 



6 F 2 -Model 

We consider the special model with Pj 



0, 1 in some detail. The boundary matrix is given by 12 



independent parameters with J2k,i = 0- The initial condition Hvq = —E2V0 is explicitly solved 



as 



yo(o) 



y (i) = 



0-00 
#00 

#10 
^0°0° 

0-00 

"01 

rrOO 
#10 



trlO 
#00 
rrlO 
#10 



#00 

H\\ + 1 
- 1 



H° l + 1 
0-01 

-"10 



trlO 

-"oo 

trlO 
-"01 
rrlO 
-"10 



(45) 



(46) 



Upon setting the function m(t) as 



m(t) 



fj-OO / 
-"00 - * 


0-01 

#00 


rrlO 
#00 


#00 


yo(o) 


#o°i° + 1 


#0°l 


rrlO 
-"01 


-"01 


-Vo(0) 


rrOO 
#10 


trOl 
#10 


#10 - t 


#10 


yo(i) 


-t 


1-t 








yo(o) 








-t 


1 - t 


yo(i) 



(47) 



the condition is shown m(t) = independent of t, which gives 3 relations among the 12 parameters H* 3 
As an example, we concentrate the model with 



kl- 



trOl 
#11 



a, H™ = b, 



rrlO 
#01 



p, H§=q, H$ = r, 



00 



and the other parameters = for / (A;, I). The equation gives the relation 

a + q = ap, (r + b) (ap — r — b) + pr = 0. 



(48) 



For instance, the model with a = b = l,p = 2,q = 1 and r = 1 + y/2 satisfies this condition. It is obvious 
that the simplest ASEP model (a = b,p = 1 and q = r = 0) is the simplest example. 

Let us consider the eigenvalue problem. In this case, we see that there are two eigenvalues (/ = 2) as 
well as the simplest ASEP model and the eigenvalues a± and a^ are calculated as 



ai 



b + r' 



(l 2 = 1. 



Therefore, we see that if the parameters satisfy 1/(6 + r) < 2(> 2), the system indicates the maximal 
current phase (low density phase). 
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7 Summary and Discussion 



As we have seen, possible integrable boundary interaction should satisfy eq. (pl|). The density current and 
concentration can be expressed in terms of the partition functions exactly determined by the recursion 
relations The extension for general /3 and 7 at the right edge and p and q in the bulk is 

straightforward. Because stochastic systems are equivalent to ferromagnetic spin chains, their behavior 
highly depends on the boundary conditions. 

We must remark several points. First, the fact that an eigenvalue a<± equals to 1 originates in the 
choice that we take an extraction rate 1 at the right edge. Secondly, some of the readers may wonder; 

If we hide or ignore the first site, does the system reduce to the simplest A SEP model? 

In fact, the answer is no. There are two reasons or evidences. First, we must pay attention that a\ is 
not equal to just the sum of the probability YUj Hji because of the conservation of probability at the 
left edge. Secondly, even if we focus on the simplest ASEP model, we cannot obtain the precise results 
by hiding the left edge site ||, 1C]. This fact comes from the recursion relations (|7|)- (|9|). The form of a\ 



indicates an effect of bulk correlation, though it is not so clear because we have set p = 1 in this case. 

In general, the current J of the system does not depend on the parameters H^ p a explicitly but on 
the eigenvalues cij. It is still unknown whether there exists an interesting model with X G V s (s > 3), 
especially XeV^(g>V^<g>---(g)V^. 
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